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Abstract. We report the discovery of very-high-energy (VHE) y-ray emission of the binary system PSRB 1259-63 / SS 2883 
of a radio pulsar orbiting a massive, luminous Be star in a highly eccentric orbit. The observations around the 2004 periastron 
passage of the pulsar were performed with the four 13 m Cherenkov telescopes of the H.E.S.S. experiment, recently installed 
in Namibia and in full operation since December 2003. Between February and June 2004, a y-ray signal from the binary 
system was detected with a total significance above 13 cr. The flux was found to vary significantly on timescales of days 
which makes PSRB 1259-63 the first variable galactic source of VHE y-rays observed so far. Strong emission signals were 
observed in pre- and post-periastron phases with a flux minimum around periastron, followed by a gradual flux decrease in 
the months after. The measured time-averaged energy spectrum above a mean threshold energy of 380 GeV can be fitted by a 
simple power law F (E/l TeV)~ r with a photon index T = 2.7 ± 0.2 stat ± 0.2 sys and flux normalisation F = (1.3 ± 0. l sta t ± 
0.3 sys ) x 10~ 12 TeV~' cnT 2 s -1 . This detection of VHE y-rays provides unambiguous evidence for particle acceleration to multi- 
TeV energies in the binary system. In combination with coeval observations of the X-ray synchrotron emission by the RXTE 
and INTEGRAL instruments, and assuming the VHE y-ray emission to be produced by the inverse Compton mechanism, the 
magnetic field strength can be directly estimated to be of the order of 1 G. 
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1. Introduction 

PSRB1259-63/SS2883 is a binary system consisting of a 
~48 ms pulsar in orbit around a massive B2e companion star 
(Johnston et al. 1992a,b). The highly eccentric orbit of the pul- 
sar places it just ~ 10 13 cm from the companion during perias- 
tron every ~3. 4 years. Be stars are known to have non-isotropic 
stellar winds forming an equatorial disk with enhanced mass 
outflow (e.g. Waters et al. 1988). In the case of PSRB1259-63, 
timing measurements suggest that the disk is inclined with re- 
spect to the orbital plane (Wex et al. 1998), probably because 
the neutron star received a substantial birth kick, causing the 
pulsar to cross the disk two times near periastron. These unique 
properties make the binary system PSRB1259-63 an excel- 
lent laboratory for the study of pulsar winds interacting with 
a changing environment in the presence of an extremely in- 
tense photon field. The synchrotron origin of optically thin 
unpulsed radio emission detected from this source, in partic- 
ular during the periastron passage (e.g. Johnston et al. 1999; 
Connors et al. 2002), indicates acceleration of electrons to rel- 
ativistic energies. The acceleration process has been argued to 
be most efficient when the pulsar passes through the equato- 
rial disk (Ball et al. 1999). The basic features of such a sys- 
tem ("binary plerion") in the context of higher energy X- and 
y-radiafion components have been comprehensively discussed 
by Tavani & Arons (1997). 

The intense photon field provided by the companion star 
not only plays an important role in the cooling of relativistic 
electrons but also serves as perfect target for the production 
of high energy y-rays through inverse Compton (IC) scatter- 
ing (Tavani et al. 1996; Kirk et al. 1999; Ball & Kirk 2000; 
Ball & Dodd 2001; Murata et al. 2004). Some of these emis- 
sion models predict wind powered shock acceleration of elec- 
trons to multi-TeV energies, radiating predominantly through 
the synchrotron and IC channels, with the main energy release 
in the X- and high energy y-ray bands, respectively. 

The unpulsed non-thermal X-ray emission detected from 
PSRB 1259-63 throughout its orbital phase in 1992 to 1996 by 
the ROSAT and ASCA satellites (Cominsky et al. 1994; Kaspi 
et al. 1995; Hirayama et al. 1996) generally supports the syn- 
chrotron origin of X-rays. The spectrum of the synchrotron ra- 
diation seems to extend to hard X-rays/low energy y-rays as 
shown by OSSE (Grove et al. 1995) and recently confirmed by 
observations with the INTEGRAL satellite (Shaw et al. 2004). 

Since the companion star provides the dominant source of 
photons for IC scattering, the target photon density is well 
known throughout the entire orbit. Therefore, the ratio of X- 
ray flux to high energy y-ray flux depends only on the strength 
of the ambient magnetic field. Although the latter can be es- 
timated within a general magneto-hydrodynamic treatment of 
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the problem, it contains large uncertainties which affect the es- 
timate of the IC y-ray flux as F y oc B 2 . 

Kirk et al. (1999) studied the light curves of very-high- 
energy (VHE) y-rays under the assumption of a 1/r depen- 
dence of the magnetic field which implies that the ratio of the 
energy density of the photon field to that in the magnetic field 
B is independent of orbital phase. They also assumed that the 
position of the termination shock as well as the strength of the 
magnetic field is not affected by the disk of the B2e star. Under 
such assumptions, they predicted an asymmetric y-ray light 
curve with respect to periastron (because of the inclination of 
the orbit with respect to the line of sight and the dependence of 
the inverse Compton y-ray emissivity on the scattering angle), 
with an increase towards periastron and monotonic decrease af- 
ter the passage of periastron. However, one might possibly ex- 
pect significant deviation from such a simplified picture given 
the apparent strong impact of the disk on the pulsar wind termi- 
nation as seen in the X-ray light curve (Tavani & Arons 1997). 
Moreover, during the time periods of interaction of the pulsar 
wind with the equatorial disk one may expect, in addition to 
the IC y-rays, a new component of y-radiation associated with 
interactions of accelerated electrons and possibly also protons 
with the dense ambient gas (Kawachi et al. 2004). Up to now, 
the theoretical understanding of the properties of this complex 
system, involving pulsar and stellar winds interacting with each 
other, is quite limited because of the lack of constraining obser- 
vations. 

Nevertheless, the fortunate combination of: (1) the high 
spin-down luminosity of the pulsar, L = 8.3 x 10 35 erg/s, which 
is partially converted into populations of ultra-relativistic par- 
ticles, (2) the presence of the intense target photon field 
provided by the companion star with energy density « 
0.9 / 10 13 cm) ~ 2 erg/cm 3 (where R(t) is the spatial separa- 
tion between the pulsar and the companion star), and (3) the 
relatively small distance to the source (d » 1 .5 kpc) makes this 
object a very attractive candidate for VHE y-ray emission. 

Previous observations of PSRB1259-63 in VHE y-rays, 
apart from its periastron passage were performed using the 
CANGAROO I and CANGAROO II detectors, but did not re- 
sult in significant signals and provided upper limits at 13% 
of the flux from the Crab Nebula (see Kawachi et al. 2004, 
and references therein). The first significant detection in VHE 
y-rays based on preliminary analysis results was reported in 
Beilicke et al. (2004) a few days prior to the periastron pas- 
sage in February 2004 by the High Energy Stereoscopic System 
(H.E.S.S.) to allow for target of opportunity observations of 
other instruments. The analysis of the data from this initial de- 
tection and the data obtained in the subsequent H.E.S.S. obser- 
vation campaign on PSRB 1259-63 is presented in this paper. 

2. VHE y-ray Observations and Results 

2.1. Observations and Analysis 

The observations from February to June 2004 were performed 
with the High Energy Stereoscopic System (H.E.S.S.), consist- 
ing of four imaging atmospheric Cherenkov telescopes (Hinton 
2004) located in Namibia, at 23°16' S 16°30' E in 1800m 
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above sea level. Each telescope has a tesselated spherical mir- 
ror with 13 m diameter and 107 m 2 area (Bernlohr et al. 2003; 
Cornils et al. 2003) and is equipped with a camera of 960 0.16°- 
photomultiplier tubes providing a total field of view of 5° in 
diameter (Vincent et al. 2003). During the stereoscopic obser- 
vations, an array trigger requires the simultaneous detection of 
air-showers by several telescopes at the hardware level, allow- 
ing a suppression of background events (Funk et al. 2004). 

All observations were carried out in moonless nights track- 
ing sky positions with an alternating offset of typically +0.5° 
in declination relative to the source (the wobble mode) in time 
intervals of 28 minutes duration. This allows to determine the 
background from the same field of view and one can omit off- 
source observations, effectively doubling the observation time. 
Due to the serendipitous discovery of another source in the field 
of view around PSR B 1259-63 (see Aharonian et al. 2005a), 
for all observations subsequent to the 14th of May 2004 (MJD 
53139) the array pointing was changed to a position ~ 0.6° 
north of PSR B 1259 -63 and an alternating wobble offset of 
0.5° in right ascension (instead of declination) was used. 

The data set, selected on the basis of standard quality crite- 
ria, has a dead time corrected exposure (live time) of 48.6 h and 
a mean zenith angle of 42.7°. The corresponding mean thresh- 
old energy defined by the peak y-ray detection rate for a source 
with a Crab-like spectrum (Aharonian et al. 2000) after selec- 
tion cuts was estimated to be 380 GeV. In the phase prior to 
the periastron passage (7.8 h live time), due to technical prob- 
lems, data from only three telescopes were considered, and for 
the post-periastron phase (41.9 h live time) data from the full 
telescope array were used. 

After the calibration of the recorded air shower 
data (Aharonian et al. 2004), each telescope image was 
parametrised by its centre of gravity and second moments 
(Hillas 1985) followed by the stereoscopic reconstruction of 
the shower geometry providing an angular resolution of ~ 0. 1 
for individual y-rays. The y-ray energy was estimated from 
the image intensity and the shower geometry with a typical 
resolution of ~ 15%. In order to reject the vast background 
of cosmic-ray showers, y-ray candidates are selected using 
cuts on image shape scaled with their expectation values 
obtained from Monte Carlo simulations. The cuts used for 
this analysis were optimised on simulations of a y-ray point 
source with a flux level of 10% of the Crab Nebula VHE 
y-ray flux, allowing 41.2% of the y-rays to be retained while 
rejecting more than 99.9% of the cosmic-ray air showers. A 
more detailed description of the analysis techniques can be 
found in Aharonian et al. (2005b). 

2.2. Detection 

Figure 1 shows the distribution of the squared angular distance 
8 2 of excess events relative to the position of PSR B 1259-63 
for the whole data set. The circle with a radius corresponding 
to the angular cut (% ut = 0.02 deg 2 in the field of view around 
the source position was considered as the on-region. The back- 
ground was estimated from several non-overlapping circles of 
the same radius (off-regions) with the same angular distance 




0.02 0.04 0.06 0.08 0.1 

e 2 [deg 2 ] 



Fig. 1. Main figure: Distribution of background-subtracted y-ray can- 
didates in 9 2 where 8 is the angular distance between the recon- 
structed shower direction and the pulsar position for the complete 
2004 data set. The solid line indicates the distribution expected for 
a point source of y-rays. Inset: Distributions in 9 1 for "on-source" y- 
ray candidates relative to the source position (full points) and for all 
background control regions (histogram, scaled with the background 
normalisation factor a). An uniform background results in a flat distri- 
bution in 2 . The dashed vertical line indicates the angular cut applied 
(0 2 <Ct = 0.02 deg 2 ). 

from the camera centre, allowing corrections due to the vary- 
ing camera acceptance to be omitted (Aharonian et al. 2001). 
The clear excess in the direction of the pulsar has a significance 
of 13.8 cr and is consistent with a distribution obtained from a 
simulated y-ray point source. The y-ray signal from the direc- 
tion of PSR B 1259-63 was detected in most of the darkness 
periods from February to June 2004, for which the results are 
summarised in Table 1 . 

A 2D-analysis of the H.E.S.S. field of view around 
PSR B 1259-63 was performed using the ring background tech- 
nique as alternative background estimation method. In this 
method, for each bin, the number of on-events N on was de- 
rived by integrating the bin content within a circle of radius 

Table 1. Results of H.E.S.S. observations on PSR B 1259-63. For 
each darkness period the number of telescopes used in the analysis 
iV te i, the live time f H ve, the significance S calculated according to Li 
& Ma (1983), the number of counts within the on-source (N oa ) and 
off-source (AW) region(s), the background normalisation a, and the 
number of detected y-rays (N y ) are listed. The background normali- 
sation a was determined by the number of off-regions, which do not 
coincide with other y-ray sources in the field of view. Therefore a is 
not the same for all data subsets because it depends on the distribution 
of wobble offsets used when obtaining the data. 



Period 


Mel 


hive 


S 


A on 


N oS 


a 


N y 


2004 




M 


[cr] 










February 


3 


7.9 


9.8 


691 


2591 


0.172 


246±25 


March 


4 


17.8 


7.2 


1740 


8459 


0.169 


307+43 


April 


4 


5.1 


7.8 


644 


2705 


0.167 


193+25 


May 


4 


9.9 


5.2 


910 


4499 


0.167 


160+31 


June 


4 


9.1 


1.8 


722 


4027 


0.167 


51+28 


Total 




49.8 


13.8 


4707 


22281 


0.168 


955±69 



4 
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Fig. 2. Significance sky-map centred on the position of 
PSR B 1259-63 for the H.E.S.S. February data. The excess ~ 0.6° 
north of the pulsar is the newly discovered unidentified TeV source 
HESS J1303-631 (Aharonian et al. 2005a). The bins are correlated 
within a circle of radius CU( = 0.02 deg 2 . The background for each 
position in the field of view was estimated from a ring around this 
position correcting for the radial acceptance of the instrument within 
the field of view and background contamination from the two y-ray 
sources. 

#cut (on-region). The number of background events N Q s was 
estimated from a ring around the bin position with mean radius 
0.6° and an area 7 times larger than the area of the on-region. 
The normalisation a was corrected for the decrease of the ra- 
dial acceptance of the cameras towards the edge of the field of 
view. Figure 2 shows the significance sky-map of a 2.3° field of 
view around PSR B 1259-63 for the February data. The excess 
at the pulsar position (7V y = N on - aN Q tf - 227 ± 27, with a 
significance S = 9.1 cr) is consistent with the excess given in 
Table 1. The additional broad excess ~ 0.6° north of the pulsar 
is the unidentified TeV source HESS J1303-631 discovered in 
the same field of view (Aharonian et al. 2005a). The resulting 
bias in the background estimation due to both sources was cor- 
rected in the analysis by excluding events from a circle with 
radius 0.4° around each source from the background estima- 
tion. 

In order to derive the position of the pulsar excess and to 
check for possible source extension, the data was reanalysed 
using hard cuts requiring a minimal camera image intensity of 
200 photo-electrons which significantly improves the angular 
resolution and drastically reduces the cosmic ray background 
at the expense of a higher energy threshold of 750 Ge V. The 
uncorrected two dimensional excess distribution for the whole 
data set was fitted assuming a radially symmetric, Gaussian 
source intensity profile 

1(0) oc 

convolved with the point spread function of the detector. A pos- 
sible bias from the second source within the field of view was 



corrected for by including an additional non-symmetric two di- 
mensional Gaussian excess distribution in the fit, centred at the 
position of HESS J 1303-631 and allowing the orientation and 
width to vary freely. The resulting pulsar excess position of 
RA 13 h 2 m 49 s .3 + 2 s 3 stat , Dec - 63°49'53" ± 17" stat is consis- 
tent with the position of PSRB1259-63 (ARA = 9" ± 15" stat , 
ADec = 16" ± 17" stat ) within statistical errors. The system- 
atic uncertainty on the absolute pointing of the telescopes has 
been estimated in Gillessen (2003) to be ~ 20". A limit for the 
source extension <x s was found to be < 33" at 95% confidence 
level which corresponds to 0.24 pc at an assumed distance of 
1.5 kpc. Note that a possible source confusion with the vari- 
able hard X-ray source 1RXP J130159.6-635806 (Kaspi et al. 
1995), located 9' southeast from the pulsar, can be firmly ex- 
cluded. 

2.3. Energy Spectra 

The measured differential energy spectrum derived from all 
darkness periods with a significant detection of the source 
(February to May) is shown in Fig. 3. A power-law fit 

F(E) = dN/dE = F (j^ ) ) (1) 

of this spectrum yields a photon index F = 2.7±0.2 stat and Fq = 
(1.3 ±0.1 stat ) x lO^ctrrV'TeV- 1 . The integral flux above the 
mean threshold energy is F(E > 380 GeV) = (4.0 ± 0.4 stat ) x 
10~ 12 cm" 2 s _1 , equivalent to ~ 4.9% of the Crab Nebula flux 
(Aharonian et al. 2000) above this threshold. 

For each darkness period for which the significance of 
the observed y-ray excess exceeded 6cr (February, March, and 
April) a differential spectrum was derived (see Fig. 4) and the 
results of the corresponding power law fit according to Eq. (1) 
are listed in Table 2. Within statistical errors, there is no in- 
dication of time variability of the photon index. However, the 
changing flux normalisation Fq indicates source flux variabil- 
ity (see next section). Systematic errors were estimated to be 
<5T w 0.2 and 6F/F * 30%, dominated by the precision of 
the energy calibration of the instrument and variations in the 
atmospheric extinction of the Cherenkov light. 

Table 2. Parameters of the power law fit to the differential spec- 
trum for the different darkness periods of H.E.S.S. observations on 
PSR B 1259-63. For the periods May and June 2004, no spectrum 
could be derived due to insufficient statistics. Shown are the photon 
index T and the flux normalisation F (with statistical error only), the 
X 2 per number of degrees of freedom ^ 2 /ndf and the x 1 probability 
P x i for the power law fit of the spectrum, and the corresponding mean 
threshold energy E th (rounded to 10 GeV). 



Period 


r 


F 


x 1 ! 


p * 




2004 




[TeV -1 cm" 2 s" 1 ] 


ndf 




[GeV] 


February 


2.8 + 0.3 


2.1 ±0.3 x 10" 12 


1.8/4 


0.77 


370 


March 


3.0 ±0.6 


1.0 ± 0.2 x 10" 12 


1.8/4 


0.76 


420 


April 


2.8 ± 0.4 


2.0 ± 0.3 x 10" 12 


0.4/3 


0.93 


350 


Overall* 


2.7 ± 0.2 


1.3 ±0.1 x 10" 12 


2.3/5 


0.81 


380 



* includes all data from February to May 2004 
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Fig. 3. Differential energy spectrum dN/dE of y-rays from 
PSRB1259-63 using H.E.S.S. 2004 data from the periods with sig- 
nificant detection of the pulsar (February-May). The solid line shows 
the power-law fit to the spectrum (see also Table 2). 
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Fig. 4. Differential energy spectra of y-rays from PSR B 1259-63 for 
data of the periods February, March and April (from top to bottom). 
The straight line in each plot indicates the power-law fit to the partic- 
ular spectrum for which the results are listed in Table 2. Note, that the 
photon index T is constant within the statistical errors. 



2.4. Flux Variability 

The daily integral flux of y-rays above the mean threshold of 
380 GeV is shown in Fig. 5 (lower panel). Since the spectrum 
cannot be derived on a daily basis due to limited statistics, the 
integral flux F(E > E± = 380 GeV) was obtained by integrat- 
ing Eq. (1) assuming the spectral index of the time-averaged 
spectrum F = 2.7. The flux normalisation Fq was calculated 
using the measured excess of y-rays N y given by 



= F J d0f llve (0) j d£A eff (£,0)(y 



— f 

TeV/ 



with the daily live time fii ve (0) as a function of the zenith angle 
0, and the effective area A e g(E, 0) after selection cuts obtained 
from simulations. Note that this method does not depend on the 
energy reconstruction and so the complete sample of excess 
events can be used. An alternative method, using the energy 
reconstruction to calculate the flux based on the effective area 
for each event above the energy threshold E± according to 



F(> E,h) 




1 



1 



Aeff(£',0,) 



E'j>E th 

a y 



with E' as the reconstructed energy and a as the background 
normalisation, yields consistent results but introduces a bigger 
statistical error since only events above the threshold are con- 
sidered. 

The daily light curve in Fig. 5 clearly indicates a variable 
flux. This can be quantified by a fit of a constant flux to the data 
yielding a x 2 of 90.9 per 35 degrees of freedom corresponding 
to ax 2 probability of 1 .2 x 10" 7 . 

In order to investigate the flux trend in each darkness pe- 
riod, the individual light curves were fitted separately by a 
straight line and the significance for an increasing or decreas- 
ing flux is listed in Table 3. The low flux state after periastron 
is followed by a distinct rise beginning at r ~ +15 days and a 
slow decrease until t ~ +75 days where the observed excess is 
no longer significant. 



Table 3. Light curve properties of the H.E.S.S. 2004 data. The in- 
tegrated flux F(> 380 GeV) was calculated assuming a photon index 
T = 2.7 (see text). For the daily integral flux light curve of each period 
the following properties are given: the^ 2 /ndf for a fit to a constant 
flux, and the slope a and significance S a = \a\ /cr a of an increasing or 
declining flux after fitting a straight line F(?mjd) = a'MJD+^*(0) (?mjd is 
the Modified Julian Date), together with the corresponding ^ inc /ndf . 



Period 


F(> 380 GeV) 


X const / 


a [lO-'cm" 2 


^line / 


Sa 


2004 


fl0- 12 cm- 2 s" 1 ] 


ndf 


s" 1 MJD" 1 ] 


ndf 


[o-l 


February 


5.9 ± 0.6 


9.4/5 


-10±5 


4.3/4 


2.2 


March 


2.8 ± 0.4 


31.9/9 


7.1 ± 1.5 


8.9/8 


4.6 


April 


5.1 ±0.7 


4.2/5 


-4.4 ± 2.4 


0.8/4 


1.8 


May 


2.4 ± 0.5 


4.0/7 


-1.6 ± 1.2 


2.2/6 


1.3 


June 


< 2.0* 


4.3/5 


-0.5 ±3.1 


4.3/4 


0.2 



* 99% CL, calculated according to Feldman & Cousins (1998) 



The upper panel of Fig. 5 shows the light curve of the tran- 
sient unpulsed radio emission obtained from 2004 observations 
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Fig. 5. VHE y-ray and radio light curves of PSR B 1 259-63 around its periastron passage (epoch t=0, dotted vertical line). Upper panel, open 
points: Flux density of the transient unpulsed radio emission at 1.4 GHz (Johnston et al. 2005). The pulsed radio emission was eclipsed in the 
time interval between ~ t - 18 days and ~ r + 15 days, indicated by the two dashed vertical lines. Lower panel, full points: Daily integral flux 
above 380 GeV as measured by H.E.S.S. Note that the moonlight prohibits H.E.S.S. observations for ~ 10 days between each darkness period. 



of PSRB1259-63 (Johnston et al. 2005) for the same time 
range as for the TeV band. There is some correlation visible be- 
tween the radio and the TeV bands, especially a low flux state 
around periastron and a high flux state at ~ t + 20 days. The 
pulsed radio emission from the pulsar was eclipsed during the 
periastron passage (Johnston et al. 1992a), interpreted as due to 
scattering in the dense disk material and allowing the dates of 
the disk crossings to be estimated, in particular ~ t - 16 days 
and ~ t + 13 days for the 2004 periastron passage (Johnston 
et al. 2005). The observed dates of enhanced VHE y-ray flux 
roughly correspond to orbital phases a few days after the as- 
sumed disk crossings. However, the coverage of the considered 
time interval by H.E.S.S. observations is limited and does not 
allow firm conclusions about VHE y-ray flux maxima, espe- 
cially pre-periastron. 

3. Discussion and Conclusions 

The detection of VHE y-rays from PSRB1259-63 gives the 
first unambiguous and model-independent evidence of parti- 
cle acceleration to multi-TeV energies in this unique binary 
system. The (most likely) synchrotron origin of X-rays (e.g. 
Tavani & Arons 1997) also indicates the presence of ultra- 
relativistic particles, in the form of electrons, but other possi- 
ble explanations of the X-radiation, e.g. due to the bulk motion 
Comptonisation of the stellar photons by the pulsar wind with a 
moderate Lorentz factor of 10-100 (Chernyakova & Illarionov 
1999), cannot be firmly excluded. 



In contrast to the X-ray emission, the extension of the y- 
ray spectrum to several TeV necessarily implies that the parent 
particles, electrons and/or protons, are accelerated to at least 
10 TeV. The point-like feature of the detected signal constrains 
the extension of the y-ray production region to R < 0.24 pc = 
7 .4x 10 17 cm. The information contained in the time variability 
of the signal on timescales of days provides a hundred times 
more stringent upper limit on the size of the y-ray source of 
about 10 16 cm. 

The most likely scenario of particle acceleration and ra- 
diation in this system is a variety of the standard model of 
Pulsar Driven Nebulae (see e.g. Rees & Gunn 1974; Kennel 
& Coroniti 1984) which postulates that the deceleration of the 
ultra-relativistic pulsar wind (with a bulk motion Lorentz fac- 
tor of 10 6 — 10 7 ) by the pressure forces of the external medium 
enforces a termination shock. In this context electrons are ac- 
celerated to multi-TeV energies. These electrons radiate in the 
magnetic and photon fields in which they propagate and thus 
produce synchrotron and Inverse-Compton nebulae with a typ- 
ical size exceeding 0.1 pc, commonly called plerions. 

The mean energy flux contained in the VHE y-ray emis- 
sion derived from the time-averaged energy spectrum is EFe ~ 
3 x 10 _12 ergcm _2 s _1 which represents a y-ray luminosity of 
L y « 8 x 10 32 erg s _1 for an assumed distance of the binary sys- 
tem of 1.5 kpc. This corresponds to « 0.1% of the pulsar spin- 
down luminosity. Interestingly, this value matches the typical 
ratio between the pulsed and unpulsed component of X-ray ra- 
diation from isolated pulsars with associated pulsar wind nebu- 
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lae and their spin-down luminosity (see e.g. Cheng et al. 2004), 
supporting the suggestion that the energy for the observed y- 
rays is provided by the pulsar. 

However, in the case of PSR B 1259-63 both the spatial and 
temporal scales are expected to be rather different from those in 
plerions around isolated pulsars. Also, dynamical flow effects 
can play a different role since the pulsar moves through the 
stellar wind and the photon field of the massive companion star. 
Due to the high pressure of the stellar wind, the pulsar wind 
terminates very close to the pulsar and thus the electrons are 
accelerated well within the binary system. Moreover, because 
of severe adiabatic and radiative losses, the TeV electrons have 
a very short lifetime. As a result, the radiation is emitted in a 
rather compact region not far from the acceleration site, and for 
any given time the emission originates from a quite short sector 
of the pulsar trajectory. 

In the following we summarise the physical processes ex- 
pected to be at work, and put our observational results into this 
perspective. 

3.1. Characteristic Cooling Times 



wind interface 



The synchrotron loss time 



<sy 



6n 



2 4 

crjceB 2 



400e X( ! v flG 2 s, 



(2) 



of an electron with energy e depends on the local magnetic field 
strength B (where Be = B/IG) which is not known a priori. 
We shall estimate B below by combining the INTEGRAL X-ray 
data and our H.E.S.S. TeV y-ray data. In contrast, we can es- 
timate the inverse Compton cooling rate of electrons with very 
good accuracy at any given time, i.e. at any location of the pul- 
sar in its orbit (Tavani & Arons 1997; Kirk et al. 1999). 

Indeed, the inverse Compton losses are dominated by the 
scattering of electrons on the light of the companion star with 
luminosity L star = 3.3 x 10 37 erg/s and effective blackbody 
temperature T e g = 2.3 x 10 4 K. At any distance R, the black- 
body radiation decreases by a factor of (R/R+y 2 , where R+ = 
4 x 10 11 cm is the radius of the companion star. Thus at the 
distance of the pulsar to the star = R/IQ U cm, the energy 
density of the radiation is 



W T = 



4ttR 2 c 



0.87fl B 2 erg/cm 3 , 



(3) 



and correspondingly the electron cooling time is determined in 
the deep relativistic Klein-Nishina regime which applies here 
for TeV energies, given the high energies of the seed photons. 
For a Planckian seed photon spectrum, and using the formalism 
of Blumenthal & Gould (1970), one obtains 



?kn = 3.2 x 10- 



fTeV^i3 



ln(30.6erev)- 1-4 



(4) 



For PSR B1259-63, with its environment of a strong stellar 
wind and a strong radiative flow from the companion, we must 
in addition consider the adiabatic losses of the accelerated pul- 
sar particles. This has already been noted by Tavani & Arons 
(1997). Indeed, the plasma flow beyond the pulsar wind termi- 
nation shock differs strongly from the spherically symmetric 



periastron 
epoch x 




H.E.S.S. 

PSR B1 259-63 



pulsar orbit 
period = 3.4 years 
eccentricity = 0.87 



Fig. 6. Sketch of the orbit of PSR B 1 259-63 with respect to the line of 
sight (adapted from Johnston et al. 1999). The pulsar approaches the 
equatorial disk prior to periastron while it is "behind" the companion 
star and turns towards the observer before it crosses the disk after pe- 
riastron for the second time. Note that the orientation of the disk with 
respect to the orbital plane is not precisely known. The colour gradient 
bars along the orbit indicate the periods of H.E.S.S. observations and 
show the integral VHE y-ray flux using a smoothed light curve based 
on the data points from Fig. 5. 

configuration of an isolated pulsar in a static pressure environ- 
ment. Rather than being decelerated in a spherically symmetric 
fashion while compressing the post-shock magnetic field, the 
actual flow will stagnate at the sub-stellar point at a distance 
I ~ 10 12 cm and expand relativistically into the downstream di- 
rection, which roughly points away from the Be star. The sur- 
face of separation between the interior flow of the pulsar par- 
ticles and the exterior stellar wind flow is indicated as "Wind 
interface" in Fig. 6. 

The expansion of the relativistic, shocked pulsar wind gas 
will proceed with a flow velocity of magnitude |m|, roughly 
equal to the relativistic magnetosonic speed ev mS;J _, where 
Vms,-L - c V2/3, over a distance Al, where Al is a few xl. The 
factor s ranges between V2/2 and 1 , if the direction of the post- 
shock MHD flow ranges between a perpendicular and a parallel 
orientation with respect to the magnetic field direction, respec- 
tively. Assuming the pitch angle scattering mean free path of 
the accelerated pulsar wind particles to be small compared to I, 
their adiabatic loss time is then given by 



fad 



3e 
divM 



s x 350 



Al 



2x 10 12 cm' 



(5) 



with 0.71 < e < 1, depending on the details of the post-shock 
MHD flow. 

In this picture the adiabatic losses proceed about as fast as 
the synchrotron losses if Z? as 1 G (compare with Eq. (2)) - 
while both are faster than the inverse Compton losses - and 
the energy spectrum of the radiating electrons is close to their 
source spectrum produced during acceleration at the termina- 
tion shock. We shall consider neither electron Bremsstrahlung 
nor y-ray production in interactions of relativistic protons with 



8 



F. Aharonian et al.: Discovery of PSR B1259-63 in VHE y-rays around Periastron 



the ambient gas via decay of secondary 7r°-mesons (see e.g. 
Kawachi et al. 2004) here, since the accelerated pulsar wind 
particles are hydromagnetically constrained to the magnetic 
field lines of the pulsar outflow and will hardly interact individ- 
ually with the stellar wind particles, except in a thin boundary 
layer at the Wind interface far downstream. 

The missing quantity in our cooling considerations is the 
value of B. We shall use the contemporaneous INTEGRAL re- 
sults and our H.E.S.S. fluxes to estimate it as follows: Since 
there is little doubt that the detected TeV y-rays are produced 
in the Klein-Nishina regime, we may conclude that they are 
emitted by electrons of the same energy, e e ~ E y . This allows 
to set up a direct relation between y-rays and X-rays produced 
by the same electrons, 



E x >20fi G £ 2 eV keV. 



(6) 



The Compton losses with Ej e \ take place in a spatial region 
that is larger than that of synchrotron losses, and therefore adi- 
abatic losses have already occurred there. 

An experimental estimate of B is possible by comparing the 
energy fluxes of the detected X- and y-rays. The post-periastron 
data obtained in March by RXTE and INTEGRAL in the X-ray 
band from 1 keV to lOOkeV (see Shaw et al. 2004) show a 
relatively flat spectral energy distribution vF Y ~ E 2 dN/dE at 
the level of 3 x 10~ 11 erg /cm 2 s. The VHE y-ray flux for the 
same period was « 2 x 10~ 12 erg/s. This implies that in the 
TeV energy regime the synchrotron losses of electrons proceed 
an order of magnitude faster than the inverse Compton losses. 
Therefore, from the comparison of the synchrotron cooling 
time, given by Eq. (2), with the Klein-Nishina cooling time 
given by Eq. (4), and taking into account that two weeks after 
the periastron the separation between the pulsar and the com- 
panion star was approximately 1.5 x 10 13 cm, we find 



B* IE T ^G. 



Our experimentally determined fields of the order of 1 G 
agree rather well with the theoretical estimate of the field de- 
rived from the MHD treatment of the pulsar wind (Tavani & 
Arons 1997; Ball et al. 1999), and the value deduced from 
the observed unpulsed radio emission near periastron (Connors 
et al. 2002). 

3.2. Spectral and Temporal Characteristics of the TeV 
Radiation 

Let us assume that the accelerated electrons enter the y-ray pro- 
duction region with a rate Q(e) oc e~ ff °. As a consequence, in 
the case of dominant synchrotron or inverse Compton losses 
in the Thompson regime, a steeper electron spectrum is estab- 
lished in the y-ray production region with power-law index a = 
ao + 1, while dominating adiabatic (or Bremsstrahlung losses) 
will not change the initial electron spectrum, i.e. a = ao, both 
for ao > 1 . 

Let us conversely assume that the radiating electrons have 
a power law distribution dN/de oc e~ a . Then the inverse 



10" 



— dominant adiabatic losses 
dominant radiative losses 
• H.E.S.S. February data 



x-12 



o 
> 



10" 



LU 



10" 



10" 



synchrotron 



in" 13 ' L 

10 10 -i„ 



B = 0.32 G 




10 



10 



/B = 3.2 G 
10" 4 10" 2 10 u 10^ 

Energy [TeV] 



Fig. 7. Spectral energy distribution from Kirk et al. (1999) assuming 
the electron cooling to be dominated by radiative (dashed lines) or adi- 
abatic losses (solid lines) for a magnetic field strength of 0.32 G at an 
orbital phase 12 days prior to periastron. The spectra were calculated 
by Kirk et al. (1999) for different electron injection spectra with the 
power law indices a = 1 .4 and ar = 2.4, respectively. For illustra- 
tion, the predicted spectrum for adiabatic losses for B = 3.2 G is also 
shown (dotted line). Data points represent the measured energy flux 
from H.E.S.S. observations in February 2004. 



Compton scattering in the Klein-Nishina regime leads to a y- 
ray spectrum (Blumenthal & Gould 1970) 



dN 
(7) dE 



oc E 



-(<*+!) 



In^ + C 



(8) 



Using the H.E.S.S. results, the photon index of the radiating 
electrons can therefore be estimated from the observed y-ray 
spectrum dN/dE oc E~ Y ~< with T y = 2.7 ± 0.2 sta t by equating it 
with the expected Klein-Nishina spectrum of Eq. (8). A more 
detailed numerical evaluation of the IC y-ray spectrum shows 
that in our energy range the average photon spectral index can 
be approximated as T y « a + 0.5. Therefore we find a « 2.2. In 
the scenario in which adiabatic losses dominate, this is, within 
errors, compatible with an acceleration process giving ao = 2. 
If radiation losses dominate, this requires an injection spectrum 
significantly harder than ao = 2. 

The relevant photon spectral index range 1 .6 < F x < 1 .8 in 
the 1 < E x < lOkeV energy range from the 1994 ASCA ob- 
servations (Hirayama et al. 1996), although not coeval, would 
also indicate a rather hard spectrum of the radiating electrons 
2.2 < a < 2.6 which would again be difficult to reconcile with 
dominant synchrotron cooling. This 1994 synchrotron photon 
spectral index is consistent with a preliminary spectral fit with 
T x = 1.82 ± 0.34 derived from the coeval INTEGRAL observa- 
tions (Shaw 2004). In this paper we shall ultimately rely on the 
H.E.S.S. spectrum. 
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Fig. 8. Comparison of the light curves predicted by different models of VHE y-ray emission with H.E.S.S. data. The solid line was taken 
from Kawachi et al. (2004), for a misaligned disk geometry, model (iii) in Fig. 7, and scaled by a factor of 0.3. The dashed and dotted lines 
represents the curves given in Kirk et al. (1999), for dominant adiabatic and radiative energy losses, respectively. The full points represent the 
corresponding light curve obtained from the H.E.S.S. data using a bin width of 2 days. 



Within the scheme of dominant adiabatic cooling, the 
strong variability of the VHE y-ray emission might be explain- 
able by a varying spatial confinement of the accelerated pulsar 
wind particles by the kinetic and thermal pressure of the stellar 
mass outflow. In Figure 6 the variations of the measured VHE 
y-ray flux are illustrated within the context of the orbital pa- 
rameters and environment of the binary system with respect to 
the line of sight. The enhanced mass outflow in the equatorial 
disk will result in a more compact emission region - and there- 
fore weaker adiabatic cooling - during the pulsar disk passages 
compared to the other regions of the orbit, in particular the pe- 
riastron phase, where the expected stronger adiabatic cooling 
could lead to a minimum of the y-ray flux. A second flow aspect 
would be a Doppler modulation due to the relativistic bulk flow 
velocity away from the stagnation point of the shocked pulsar 
wind. Considering Fig. 6, this could again lead to a flux min- 
imum during periastron. However, a detailed consideration of 
this overall physical picture is beyond the scope of the present 
paper. 

A minimum in the X-ray flux, observed for previous peri- 
astron passages, has been interpreted as being due to increased 
IC losses (Tavani & Arons 1997). However, H.E.S.S. measure- 
ments indicate a minimum in the IC flux around periastron, 
which would rule out such an interpretation for this periastron 
passage. Unfortunately, due to the full moon and bad weather 
conditions we were not able to take data during the periastron 
passage of the pulsar. 

The predicted pre-periastron spectrum from Kirk et al. 
(1999) is compared to the corresponding H.E.S.S. data in 
Fig. 7. The model curves, taken from Kirk et al. (1999) 
(Figs. 2,6), were computed for the two scenarios of dominant 
radiative losses and dominant adiabatic losses, and matched the 
archival X-ray data for electron injection spectra with ao = 1 .4 
and ffo = 2.4, respectively. They agree well with the TeV data, 
assuming a magnetic field strength of B = 0.32 G supporting 
the estimate of the order of 1 G obtained from the H.E.S.S. 
data. Despite the spectral agreement, the models fail to describe 



the TeV light curve as can be seen in Fig. 8 (dashed and dot- 
ted lines). In particular, neither the minimum around nor the 
high flux states before and after periastron are described by the 
model. This is perhaps not surprising, because the model as- 
sumed that the ratio between the energy densities of the mag- 
netic field and the photon field does not change throughout 
the full orbit, and that the effects related to the pulsar passage 
through the stellar disk can be neglected. 

Interestingly, the observed light curve seems to be qualita- 
tively similar to the prediction made by the model of Kawachi 
et al. (2004) which is also shown in Fig. 8 (solid line). In 
this model, hadronic interactions and 7r° production in the mis- 
aligned stellar disk plays a dominant role in the y-ray produc- 
tion mechanism. However, present data do not allow safe con- 
clusions concerning the interpretation of the TeV light curve. 
Studies of the most interesting part of the pulsar orbit have to 
be postponed until multi-wavelength observations during peri- 
astron become possible. 

In summary, the detection of VHE y-ray emission from 
the binary pulsar PSR B 1259-63 by H.E.S.S. provides the first 
model-independent evidence of particle acceleration in this ob- 
ject. The results clearly demonstrate the power of y-ray obser- 
vations for the study of the properties and the nature of high 
energy processes in this unique cosmic accelerator. Further ob- 
servations at VHE y-ray energies are necessary to derive spec- 
tra on a daily basis also during periastron passage. 
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